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This paper reports 139 La NMR measurements of a powder sample of the colossal magnetoresistance 
compound La 2 /3Ca 1 / 3 Mn03 (T c =268K) performed in the paramagnetic state (292K-575K) and 
in high magnetic fields (2.00-9.40 Tesla). Analysis of the spectrum measured at 575K establishes 
that the spectrum is a standard powder pattern broadened to a significant degree by a variation in 
lattice distortions around lanthanum nuclear sites. At lower temperatures, but still above T c , the 
spectrum shifts and broadens. Both the shift and broadening exhibit Curie- Weiss behavior, indicat- 
ing that the shift measures the polarization of the electron spin system, and the broadening reflects 
a distribution of magnetic susceptibilities. This distribution may result from variations of local 
susceptibility in the bulk of the sample or from differences in demagnetizing factors among powder 
grains. Close inspection of the spectrum indicates that the lattice distortions do not change as the 
temperature lowers. Spectral diffusion measurements suggest that the temperature dependence of 
the spectrum shape does not result from freezing out of motion of magnetic polarons. Variations 
in the nuclear spin-lattice relaxation across the spectrum indicate that magnetic fluctuations, not 
lattice vibrations, dominate nuclear relaxation. Nuclear spin-lattice relaxation therefore measures 
electron spin dynamics in this system. The magnetic field dependence of the spin-lattice relaxation 
indicates that the electron spin-spin correlation function adopts simple single exponential behavior 
with a slow field-independent correlation time of 10 -8 seconds near T c . The spin-spin correlation 
function changes form at higher temperatures and can be described by introducing a field depen- 
dence to the correlation time and to the magnitude of the fluctuating field. Even at the highest 
temperatures, the correlation time remains slow, on the order of 10 -9 seconds. The spin-lattice 
relaxation therefore indicates the presence of extremely slow dynamics above T c . 



I. INTRODUCTION 

Colossal magnetoresistance manganites provide an op- 
portunity to study systems in which strong coupling 
among charge, spin, orbital, and lattice degrees of free- 
dom determine bulk properties jjj . Many experiments in- 
dicate the presence of local entities in these compounds. 
Neutron scattering experiments find small magnetic clus- 
ters in the paramagnetic state for which the correlation 
length does not diverge at T c p]. A large diffusive peak 
seen by inelastic neutron scattering [||J|] , and a broad dis- 
tribution of muon spin relaxation rates || indicate that 
clusters persist far below the Curie temperature. The 
dynamics and trapping of these clusters may be respon- 
sible for the large magnetoresistance and other unusual 
transport and magnetic properties. Neutron pair distri- 
bution function [f7 L p [ and X-ray absorption fine structure 
measurements PJ10| demonstrate the presence of local 
distortions of MnOg octahedra in which the Mn-0 bond 
lengths differ significantly from the average value. While 
these distortions change during the nearly simultaneous 
ferromagnetic and metal-insulator transitions, they per- 
sist into the metallic regime. The precise role these clus- 
ters and distortions play in the behavior of the mangan- 
ites remains unclear. 

Because characterization of local properties is needed 
to understand the manganites, information from a lo- 



cal probe such as Nuclear Magnetic Resonance (NMR) 
should be valuable. NMR measurements determine prop- 
erties at a well-defined point in the system- the nu- 
clear site- but do not require that the system exhibit 
long range order. In fact, NMR can often character- 
ize the degree of disorder in a system. Analysis of the 
NMR spectrum determines the size and distribution of 
static local magnetic susceptibilities and lattice distor- 
tions. NMR spin-lattice relaxation rate measurements 
characterize extremely low energy fluctuations (on the 
order of O.l^eV) of electron spins or of the lattice. 

This paper presents 139 La NMR spectrum and nuclear 
spin-lattice relaxation measurements of 
La2/3Ca!/ 3 Mn03 (T C =268K). Measurements were per- 
formed in the paramagnetic state (292K-575K) in high 
magnetic fields (2.00-9.40 Tesla). The sample is a powder 
synthesized by a standard solid state procedure. Chem- 
ical analysis found no traces of contaminants. T c was 
determined by magnetization measurements. The mag- 
netization transition is sharp, the full width at half max- 
imum of the ^ curve at T c in a 200 Oe field being 8K. 
The width of the transition compares favorably with that 
measured in other ceramic samples and indicates that 
the sample is single phase with small deviation from the 
main composition. 

The paper begins with an introduction to general prop- 
erties of a 139 La NMR spectrum. An analysis of the high 
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temperature data demonstrates that the spectrum is a 
standard powder pattern, broadened to a significant de- 
gree by a variation in lattice distortions around different 
lanthanum sites. As the temperature lowers towards T c , 
the spectrum shifts and broadens due to the onset of 
ferromagnetism. The shift is proportional to the electron 
spin polarization while the broadening reflects a distribu- 
tion of susceptibilities. Close inspection of the spectrum 
indicates that the lattice distortions do not change in the 
temperature range studied. Spectral diffusion measure- 
ments show no signs of the freezing out of motion of mag- 
netic polarons. After an introduction to some of the prin- 
ciples of nuclear spin relaxation, we provide evidence that 
magnetic fluctuations, and hence electron spin dynamics, 
cause the relaxation. Finally, the field dependence of the 
nuclear spin-lattice relaxation reveals unusual and slow 
electron spin dynamics in the paramagnetic state. 

II. MEASUREMENTS AND DISCUSSION 

The 139 La nucleus has a spin 1=7/2 and a mag- 
netic dipole moment fi = TnW, where j n is 
the gyromagnetic ratio of the nucleus. For 139 La, 
139 7/27r=6.014MHz/Tesla. The interaction between the 
magnetic dipole moment and a magnetic field, B, creates 
21+1 energy levels with energies E m = m/y n hB, where m 
may adopt the values I, I — 1, / — 2, ... , —I. NMR mea- 
surements detect transitions between energy levels m, m' 
such that \m — m'\~l with frequencies, loq, proportional 
to the magnetic field at the nucleus: 

^0 = JnB. (1) 

Because the gyromagnetic ratio is known to a high de- 
gree of accuracy, the resonant frequency measures the 
magnetic field at the nuclear site. 

In a paramagnet, the field, B, at the nucleus is given 
by B = (1 + K)Ho, where K, the fractional frequency 
shift of the spectrum, is proportional to the magnetic 
susceptibility at the nuclear site, and Hq is the applied 
magnetic field. Measurement of the shift of the spectrum 
thus yields, within a constant of proportionality, the lo- 
cal magnetic susceptibility. The precise relation between 
the shift and the susceptibility depends on the properties 
of the system. The shift may result from intrinsic prop- 
erties, such as the electron hyperfine field, and from bulk 
effects, such as demagnetizing fields. A distribution of 
shifts broadens the spectrum by an amount proportional 
to the product of the width of the distribution and the 
magnetization. The distribution may result from a vari- 
ation in electronic structure among sites, variations in 
demagnetizing fields among differently shaped grains in 
a powdered sample, or, in a randomly oriented powder, 
from anisotropy of the shift. 

The 139 La nucleus also possesses an electric quadrupole 
moment, Q=20 x 10 _30 meter 2 , which interacts with elec- 
tric field gradients. Electric field gradients result from 



the distribution of charges around the nucleus. The elec- 
tric quadrupole interaction creates structure in the NMR 
spectrum. A site with cubic symmetry has an electric 
field gradient of zero, so the size of the effect of the elec- 
tric quadrupole interaction on the NMR spectrum mea- 
sures the degree of deviation from cubic symmetry of the 
surroundings of the nuclear site. 

For the measurements described here, the electric 
quadrupole interaction may be considered a perturba- 
tion of the magnetic interaction. Under the influence 
of the magnetic and the electric quadrupole interaction, 
the NMR spectrum exhibits several sharp peaks in an 
ordered single crystal. The frequency of the peak corre- 
sponding to the transition between level m and m-1 is: 

3 e 2 aO 

u m = 7n(l + K)H + - 4hI ^ * (3 cos 2 9 - l)(2m - 1). 

(2) 

The first term results from the magnetic interaction. The 
second is a shift due to the electric quadrupole interac- 
tion taken to first order in perturbation theory, e is the 
charge of the electron, q is the strength of the electric 
field gradient. 9 is the angle between the principal axis 
of the electric field gradient and the magnetic field. This 
formula applies if the electric field gradient is axially sym- 
metric. 

According to the above equation, the transition be- 
tween the 77i — 1/2 and m — —1/2 levels is unaffected by 
the electric quadrupole interaction to first order. This 
transition is, however, shifted to second order by an 
amount given by: 

( 2 ) _ 9 4/ 2 f e 2 gQ \ 2 (1 - 9 cos 2 fl)(l - cos 2 9) 
Aw i/2" 64 2/ - 1 \ fi ) ln (l + K)H 

(3) 

The electric quadrupole interaction affects the transi- 
tion between the m — 1/2 and m = —1/2 levels only 
slightly but shifts the other transitions by larger amounts 
to higher and lower frequencies. The +1/2 to —1/2 tran- 
sition is known as the "central transition" while the other 
transitions are known as the "satellites." The magnetic 
field dependence of the quadrupole interaction is partic- 
ularly important. Analysis of the spectrum in this paper 
utilizes the fact that the first order shift is independent 
of magnetic field whereas the second order shift varies 
inversely with magnetic field. 

The crystal axes determine the orientation of the elec- 
tric field gradient. In a powder sample, the crystal axes, 
and hence the electric field gradient, adopt a distribu- 
tion of angles with respect to the applied magnetic field. 
This angular distribution broadens the sharp peaks of a 
single crystal into a characteristic powder pattern which 
has been described and calculated in the literature [p+| . 
In the powder pattern, the satellites spread out consider- 
ably and overlap. The first order quadrupole interaction 
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determines the satellites' width, which is therefore pro- 
portional to the size of the electric field gradient and 
independent of magnetic field. The average of the first 
order interaction over all angles is zero, so the average 
shift of the satellites due to the first order quadrupole 
interaction is zero. The satellites do, however, exhibit a 
shift due to the magnetic interaction. 

In contrast, the central transition, feeling only a sec- 
ond order interaction, is much sharper than the satel- 
lites. The second order interaction varies as q 2 /Ho, so 
the powder-average width of the central transition is in- 
versely proportional to the magnetic field. The average 
of the second order interaction over all angles is nonzero, 
so the average position of the central transition is shifted, 
also by an amount inversely proportional to the magnetic 
field. The average frequency of the central transition for 
a spin-7/2 is: 

-X/ 2 =7n(l + ^0-4(^) Mi ' k)Hq - 

(4) 

The frequency of the central transition thus contains a 
component proportional to and a component inversely 
proportional to the magnetic field. The effects of a de- 
viation from axial symmetry (a non-zero asymmetry pa- 
rameter) on the average frequency are small, being no 
more than two percent of that due to the axially sym- 
metric part of the quadrupole interaction. 

A. The Spectrum at High Temperature 

Figure |l] depicts the 139 La spectra in La 2 /3Ca 1 / 3 Mn03 
at 575K measured at 8.19 and 4.00 Tesla. The frequency 
shift is measured with respect to the unshifted resonant 
frequency derived using the value of Hq determined from 
the resonant frequency of deuterium in deuterated wa- 
ter. The spectrum consists of a narrow feature (300 kHz 
FWHM at 8.19 Tesla, 680kHz at 4.00 Tesla) and a broad 
feature (4 MHz FWHM at both fields). The inset of fig- 
ure [l] compares the narrow feature measured at 4.00 Tesla 
and 8.19 Tesla after multiplying the frequency shift by 
the value of the field. After this scaling of the frequency 
axis, the narrow features lie on top of each other, demon- 
strating that the width scales inversely with field. We can 
therefore conclude that the narrow feature corresponds 
to the central transition of the spectrum of a spin-7/2 
nucleus. Resolution of the broad feature at the lower 
field is limited by a small signal to noise ratio. However, 
the broad feature does not become wider upon lowering 
the field, indicating that it corresponds to the satellites. 
Further evidence for this assignment of the broad feature 
comes from a fit of the spectrum. 

The average shift of the central transition has a com- 
ponent which varies as q 2 /Hq. Therefore, measurements 
of the frequency of the center of the narrow feature of 



the spectrum as a function of applied magnetic field, Hq, 
determine the average strength of the electric field gra- 
dient, q. Measurements at fields ranging from 2.00 to 
8.19 Tesla determine a value of electric quadrupole cou- 
pling, vq = e ^ =26 MHz. For the sake of comparison, 
for 139 La in the layered cuprate system LaCuC>4 +( 5, vq = 
89 MHz Unlike the cuprate, the crystal structure of 
the manganite is nearly cubic, with a slight orthorhombic 
distortion [^(j. However, the nonzero electric field gradi- 
ent demonstrates the sensitivity of NMR to the deviation 
from cubic symmetry, which can be safely ignored in a 
number of other types of measurements. 

Given the measurement of the average strength of the 
electric field gradient, we can fit the spectrum. Figures ||, 
||, and U show several fits, each assuming an electric field 
gradient with axial symmetry, with insets showing ex- 
panded views of the narrow feature and of a portion of 
the wide feature. Figure || depicts a calculated powder 
pattern spectrum which uses no adjustable parameters. 
The amplitude has been chosen so that the total area 
under the calculated powder pattern is the same as the 
area under the measured spectrum. The narrow feature 
in the calculated powder pattern is the central transition, 
and the broader features are the satellites. The widths of 
the features of the calculated powder pattern and of the 
measured spectrum match closely We can also compare 
the areas under the different parts of the spectrum. In 
the calculated powder pattern, the ratio of the area un- 
der the satellites to the area under the central transition 
is 4.96. For comparison, we fit the measured spectrum to 
a superposition of a broad and a narrow gaussian. The 
ratio of the areas of the gaussians is 5.7 ± 0.5. Given the 
difficulties involved with precisely measuring absolute in- 
tensities over the large frequency range involved in this 
experiment, the area ratios show reasonable agreement. 

The theoretical spectrum shows sharp peaks which are 
not found experimentally. We can improve the fit by in- 
troducing broadening mechanisms. Figure [?] shows the 
effect of broadening by a distribution of magnetic shifts. 
The central transition of the calculated powder pattern 
matches the narrow feature of the measured spectrum 
well, but the smeared singularities of the calculated satel- 
lites (figure ||, inset on right) do not appear in the mea- 
sured spectrum. Increasing the broadening to make the 
satellites match the broad feature makes the central tran- 
sition too wide. A distribution of magnetic shifts does not 
by itself account for the spectrum shape. 

By including a distribution of electric field gradient as 
well as magnetic broadening, we can fit both the broad 
and narrow features. Figure ^ shows a powder pattern 
with a gaussian distribution of the electric quadrupole 
coupling of width Avq / vq=15%. This substantial dis- 
tribution indicates that distortions in the surroundings 
of different lanthanum sites vary significantly. The fit 
also involves a convolution with a gaussian with a width 
0.23% of the average NMR frequency which represents a 
slight distribution of magnetic shifts. A different distri- 
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bution or a deviation of the electric field gradient from 
axial symmetry would improve the fit. However, the 
quality of this simple fit confirms the interpretation of 
the spectrum as a central transition and satellites. 

The presence of a distribution in lattice distortions 
agrees with results of other local measurements of the lat- 
tice j7|-^0[ which find a distribution of manganese-oxygen 
and lanthanum-oxygen distances. To affect the electric 
field gradient at the lanthanum site, the distortions need 
not be in bonds involving the lanthanum nucleus. Varia- 
tions in distant manganese-oxygen bonds can have a mea- 
surable effect. Such distortions may result from defects, 
Jahn- Teller distortions of MnC>6 octahedra, or from dis- 
tortions due to the difference between the radii of La 3+ 
and Ca 2+ . We do not present a quantitative comparison 
of distortions seen by NMR and those measured by other 
methods because accurate calculations of electric field 
gradients require characterization of the charge distribu- 
tion around the nucleus to a level of precision beyond the 
scope of this investigation. 

One may also reasonably expect that local lattice dis- 
tortions would cause deviations from axial symmetry. 
However, including a distribution of such deviations in 
the analysis has a negligible effect on the values of the size 
and distribution of electric quadrupole coupling. This in- 
sensitivity also precludes quantitative characterization of 
the asymmetry parameter. 

Our analysis of the spectrum rules out an interpreta- 
tion of the structure of the spectrum which invokes a 
number of sharp features corresponding to distinct sites. 
In such an interpretation, the narrow feature would, for 
example, correspond to a single site of high symmetry 
while the broad feature would correspond to a site of 
lower symmetry or to a number of narrower, overlapping 
components. Such an interpretation explains the struc- 
ture of the spectrum in the ferromagnetic state Jl3| , p^ |, 
but does not apply in the paramagnetic state. 

B. Temperature Dependence of the Spectrum 

Figure || depicts the results at several temperatures 
ranging from 292K to 575K. The spectrum experiences 
a change in shape and an enormous frequency shift on 
the order of 20% at room temperature. The shift is com- 
parable to the zero field resonant frequency of 20MHz 
observed in the ferromagnetic state |13| , |14| . In contrast, 
the bulk susceptibility creates a macroscopic demagne- 
tizing field on the order of 1% which alone would lead to 
a small negative shift. The large shift and sizable zero 
field frequency arise from the electron-nucleus hyperfine 
coupling, 

H hf = I-A-S, (5) 

where / and S are the nuclear and electron spin operators 
and A is the hyperfine coupling constant. Assuming an 
isotropic hyperfine coupling, the shift adopts the form: 



where Xe is the static electron susceptibility, and j n , "f e 
are the gyromagnetic ratios of the nucleus and of the 
electron. The hyperfine coupling constant, A, provides a 
measurement of the electron wave function at the nuclear 
site. Given the hyperfine coupling constant, the shift 
gives an independent measurement of the susceptibility. 

A plot of the inverse of the shift of the spectrum center 
of mass versus temperature combined with its magnitude 
and sign confirms that hyperfine coupling to the electron 
system determines the shift. Figure || shows that the shift 
obeys Curie- Weiss behavior with a Curie temperature of 
269±19K, in agreement with the value of T C =268K de- 
termined from magnetization measurements. The large 
error bars on the Curie temperature result from the small 
number of data points and also reflect the deviation of 
this system from simple Curie- Weiss behavior H , a hall- 
mark of short range ferromagnetic correlations. Never- 
theless, the shift clearly depends on the polarization of 
the electron system. 

Along with the shift, the spectrum changes shape. As 
figure |5| shows, the narrow feature disappears as the tem- 
perature drops. However, superimposing the room tem- 
perature spectrum and the broad part of the high tem- 
perature spectrum (figure |?]) indicates that the overall 
width remains unchanged. As explained earlier, the size 
of the electric field gradient, a result of static lattice dis- 
tortions, determines the overall width of the high temper- 
ature spectrum. The similarity between the high temper- 
ature and room temperature spectra indicates that, from 
the viewpoint of the lanthanum site, the size and distribu- 
tion of lattice distortions change little with temperature 
in the paramagnetic state. 

As we explain below, the loss of the narrow central 
feature results from a distribution of susceptibility re- 
sulting either from a variation of local electronic struc- 
ture among nuclear sites or from a distribution of demag- 
netizing fields among differently shaped powder grains. 
The resulting broadening of the narrow feature is pro- 
portional to the product of the magnetization and the 
width of the distribution. To characterize the size of 
the broadening, we start from the fit of the spectrum 
at high temperatures depicted in figure The fit in- 
cludes a slight magnetic broadening calculated by con- 
volution with a gaussian. We can fit the spectrum at 
all temperatures by varying the width of the gaussian. 
A magnetic distribution sufficient to broaden the narrow 
feature to match the width of the wide feature affects the 
wide feature only slightly. The width of the convolution 
gaussian exhibits Curie- Weiss behavior (figure ||) with a 
Curie temperature of 280±30K, in agreement with the 
value found from the temperature behavior of the shift. 
Therefore, the broadening is proportional to the magneti- 
zation, consistent with a picture in which a distribution 
in susceptibility drives the temperature dependence of 
the width of the central transition. 
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Another possible explanation for the change in the 
spectrum shape involves changes in the hopping rate of 
small polarons or other local field inhomogeneities. A 
number of experiments detect the presence of small mag- 
netic polarons in the manganites [^| JT5|JT^ |. Associated 
with the small magnetic polarons are inhomogeneous 
magnetic fields which could affect the spectrum shape. 
Nuclei near the localized moment of the small polaron 
would experience a larger shift than those far away. If the 
polaron is fixed in space, its inhomogeneous field would 
broaden the spectrum. If the moment hops at a rate ex- 
ceeding the root mean squared frequency shift induced by 
its inhomogeneous field, the nuclear system would exhibit 
a single, average shift instead of a distribution of shifts. 
As a result, the spectrum would narrow with increasing 
hopping rate. This phenomenon is known as motional 
narrowing [fl7|| . 

To see if changes in the hopping rates of local mo- 
ments affect the spectrum, we look for spectral diffusion 
by performing a type of hole burning experiment called 
an S-Wave measurement, which is described in more de- 
tail elsewhere [ |l8[ . These measurements were performed 
at the peak of the central transition in a 4 Tesla field 
at room temperature, 350K, and 575K. To perform the 
measurement, we use pulsed NMR techniques to selec- 
tively invert one section of the spectrum [figure ^] and 
watch the spectrum as it evolves over a time, T e „, af- 
ter the inversion. For the sake of explanation, assume 
that the signal at the higher frequency corresponds to 
nuclei near a local moment and the signal at the lower 
frequency corresponds to nuclei far away. We start the 
experiment by inverting the signal, and hence the nu- 
clear magnetization, on the high frequency side of the 
spectrum. If the local moment hops while the spectrum 
recovers, the nuclei originally near the moment will then 
be far from the moment and hence their signal will move 
to the lower frequency. As the inverted magnetization 
also moves to the lower frequency, a dip in the spectrum 
would simultaneously appear at the lower frequency. The 
situation involving a distribution of local fields is more 
complicated, but a dip in the signal as time progresses 
generally indicates the presence of hopping. In the ab- 
sence of hopping, the entire spectrum only grows with 
time via spin-lattice relaxation. Figure ^| shows that the 
magnetization only grows over time, suggesting that hop- 
ping does not affect the spectrum at room temperature. 
Similar measurements suggest that hopping does not af- 
fect the spectrum at 350K or 575K. 

The lack of spectral diffusion in the S-Wave measure- 
ment may, however, stem from limitations in the exper- 
iment. First, the excitation bandwidth is at most one- 
fourth the width of the central transition, so effects of 
hopping on parts of the spectrum outside the excitation 
bandwidth would not be detected. Also, if there is no 
correspondence between spectrum position and the spa- 
tial separation between a nucleus and a polaron, hopping 
of a polaron in space will not lead to spectral diffusion. 
Our measurements thus do not conclusively rule out hop- 



ping as a cause for the narrowing of the spectrum with 
increasing temperature. 

Furthermore, even if we accept that small polarons do 
not contribute to the change in the spectrum shape, such 
a conclusion would not contradict other measurements 
which detect the presence of small polarons. The density 
of charge carriers, and hence of small polarons, is rather 
large compared to that found in conventional small po- 
laron systems j[s|. Furthermore, DeTeresa et. al. || find 
that an applied magnetic field comparable to those used 
in these experiments increases the correlation length of 
magnetic clusters. Due to their density and size, small 
polarons may not present an inhomogeneous field to the 
lanthanum sites. Also, the motion of small polarons may 
be too fast, over the entire temperature regime studied, 
to affect the spectrum shape. The hopping rate in the 
temperature range studied is, at the slowest, on the or- 
der of 10 11 Hz JTq] . Motional narrowing of a spectrum 
occurs when the hopping rate reaches and exceeds the 
width of the spectrum which, in this case, is on the or- 
der of 4 MHz. We should therefore expect that changes 
in the dynamics of small polarons will have no effect on 
the spectrum shape. However, the 4 MHz lower bound on 
the rate of motion of a field inhomogeneity applies to any 
type of field inhomogeneity — an assertion independent 
of models which invoke the presence of small polarons. 

C. Relaxation Times 

The NMR spectral shape yields information about 
static, or nearly static, structural properties of the sys- 
tem. Measurements of nuclear spin-lattice relaxation de- 
termine the intensity and correlation times of magnetic 
and lattice fluctuations. In typical magnetic systems, 
magnetic fluctuations dominate the relaxation. The 
manganites, however, exhibit strong dynamic lattice dis- 
tortions |2(| which induce fluctuations in the electric field 
gradient that may dominate the relaxation. As explained 
below, measurements of nuclear spin-lattice relaxation 
performed at different points along the spectrum indi- 
cate that magnetic fluctuations dominate the relaxation. 

We start with a discussion of general aspects of nuclear 
spin-lattice relaxation of a spin-7/2 nucleus. At equi- 
librium, the population of nuclear energy levels obeys 
a Boltzmann distribution. After the populations are dis- 
turbed from equilibrium, they approach their equilibrium 
values in a characteristic time known as the spin-lattice 
relaxation time. The size of the NMR signal depends on 
the population difference between levels, so the rate at 
which the NMR signal approaches its equilibrium value 
after a perturbation depends on the spin-lattice relax- 
ation time. One type of measurement of spin-lattice re- 
laxation time is the so-called inversion-recovery experi- 
ment, in which the populations of two levels are inverted. 
The signal, M(t), associated with these levels is measured 
at several times t after the inversion. The rate at which 
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M(t) approaches its equilibrium value depends on the 
spin lattice relaxation time. 

The intensity of the transition between levels m and 
m+1 of a spin-7/2 nucleus exhibits multiple- exponential 
spin lattice relaxation pl] , p2| : 

M m (Ooc^l-a™expf-^j. (7) 

The spin-lattice relaxation time, r, depends on the in- 
tensity of fluctuations in the magnetic field and electric 
field gradient at the nucleus. The constants a™, b™, de- 
pend on the initial conditions of the populations of the 
different energy levels. Depending on whether fluctua- 
tions of the magnetic field or of the electric field gradient 
dominate the relaxation, these constants vary with m dif- 
ferently. Calculations described in the appendix indicate 
that, for a spin-7/2, the multiple-exponential spin lattice 
relaxation differs only slightly from a single exponential, 

M m (i)cxl-2exp(-^). (8) 

Furthermore, as explained in the appendix, if magnetic 
fluctuations dominate relaxation, T™ is smaller when 
measured on the central transition than on the satellites. 
If electric field gradient fluctuations dominate, the oppo- 
site is true. 

Measurements of spin-lattice relaxation at different 
points along a spectrum therefore determine the relax- 
ation mechanism. Figure [ITj depicts effective single- 
exponential relaxation times measured at different posi- 
tions on the spectrum at 9.4 Tesla and at approximately 
575K. The repetition time of this and all subsequent mea- 
surements has been set sufficiently long to ensure that 
the spin system starts at equilibrium before each mea- 
surement. The spectrum at 9.4 Tesla is similar to that 
shown in Figure [j], except that the width of central fea- 
ture (265 kHz FWHM) is smaller at higher field. The 
values of relaxation times have been normalized to the 
value measured at the apex of the central transition, and 
the frequency shift refers to the frequency of the apex. 
The spin-lattice relaxation time is shortest at the apex 
of the central transition, indicating that magnetic fluctu- 
ations cause the relaxation. 

Quantitative analysis of the variation of the relaxation 
at different points along the spectrum depends on details 
of the spectrum shape. Our analysis of the spectrum indi- 
cates that the various transitions overlap. Measurements 
of spin-lattice relaxation at any particular frequency will 
therefore contain contributions from several transitions 
at once. From the parameters used in the fit to the spec- 
trum shown in figure ^, we can determine the degree 
of overlap among transitions and so can calculate how 
the spin-lattice relaxation time varies among the differ- 
ent transitions. There is an average spin-lattice relax- 
ation time due to contributions from overlapping transi- 
tions. The line in figure [TT^ represents the calculated aver- 
age spin-lattice relaxation time as a function of position 



along the spectrum. The agreement between calculation 
and measurement supports our interpretation of the ori- 
gin of spin-lattice relaxation. Since this calculation relies 
on the parameters used in the fit of the spectrum shape, 
the agreement also lends further support to our analysis 
of the spectrum shape. 

Nuclear spin-lattice relaxation measurements probe 
magnetic fluctuations and hence electron spin dynamics. 
The spin-lattice relaxation time is related to the electron 
spin-spin correlation function via: 

- oc J dtcosuj t < h(t)h(0) > (9) 

where u>o is the NMR frequency and h(t) is the fluctu- 
ating field at the nucleus due to the electron spin. The 
constant of proportionality depends on the magnetic mo- 
ment of the nucleus and on the transition excited in the 
experiment. One consequence of the above formula is 
that the nuclear spin-lattice relaxation time depends on 
fluctuations at the NMR frequency, which are of extraor- 
dinarily low energy, on the order of 0.1/xeV in this system. 
For the sake of comparison, neutron scattering experi- 
ments can resolve fluctuations down to the meV range. 
NMR as a probe of spin dynamics thus complements in- 
elastic neutron scattering. 

If we assume a simple exponential form for the electron 
spin-spin correlation function: 

< h(t)h(0) >= hi exp (-^r) > ( 10 ) 

in which ho is the magnitude of the fluctuating field and 
t c is the correlation time of the fluctuation, the spin- 
lattice relaxation time adopts the form [^3| : 

1 2u2 u2 

1 _ 7rt"0 ll T c /-.-.x 
T 3 1+ U^c 

where the value of the constant applies for measurements 
at the central transition for a spin-7/2 nucleus. Assum- 
ing that ho and r c are field- independent, the spin-lattice 
relaxation time varies linearly with the square of the ap- 
plied magnetic field, Hg . Figure |ll| shows that, at room 
temperature, the field dependence of the spin-lattice re- 
laxation time exhibits such behavior. Furthermore, we 
find ho — 130 Gauss and r c = 10 -8 second. These values 
are only slightly affected by the presence of overlapping 
transitions in the powder pattern spectrum. The size of 
the fluctuating field is the same order of magnitude as 
the dipolar field of an electron sitting on a near neigh- 
bor manganese site, but is several orders of magnitude 
smaller than the hyperfine field which describes the shift. 
The correlation time is also several orders of magnitude 
slower than that associated with small polaron hopping. 
The slow dynamics may be related to that seen in the fer- 
romagnetic state by muon spin relaxation measurements 
Q which imply spin glass-like behavior of magnetic clus- 
ters. Electron spin resonance measurements, however, do 
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not see anomalous line broadening above T c which one 
would associate with glassy dynamics . 

At higher temperatures, the field dependence of the 
nuclear spin-lattice relaxation time changes. At 375K, 
the spin-lattice relaxation time is proportional to Hq, 
not Hq. The simplest possible explanation for such be- 
havior involves the same correlation function used above 

— 1/2 

but with field dependent parameters: r c oc H and 
h oc H^ A . Due to limits on the signal to noise ratio, 
we cannot determine whether the spin-lattice relaxation 
time at highest measured temperature, 575K, varies as 
Hq or Hq. However, regardless of the form used to fit the 
data, the correlation time at the higher temperatures is 
faster than that at room temperature, but still slow, on 
the order of 10 -9 seconds. 

In a typical ferromagnet, the magnitude of fluctuating 
fields can adopt a wide range of values. For example, 
the field from spins on the manganese may add construc- 
tively or destructively at the lanthanum site, leading to a 
random distribution of fields among different lanthanum 
sites. In the case where the size of the fluctuating field 
may adopt a range of values, the electron spin-spin corre- 
lation function should adopt the form of a gaussian. An 
exponential spin-spin correlation function results if the 
fluctuating field can adopt one of only two values. In the 
case of a gaussian correlation function, the spin-lattice 
relaxation time adopts the form: 

ln(r) K« + bH% (12) 

which does not fit any of the data, except possibly at the 
highest temperature, where the uncertainties in the data 
are large. The field dependence of nuclear spin-lattice 
relaxation thus reveals an unusual functional form for 
the electron spin-spin correlation function as well as un- 
usually slow dynamics. The explanation for this unusual 
behavior remains unclear. 



III. CONCLUSIONS 

We have established the properties of 139 La NMR of 
La 2 /3Ca 1 / 3 Mn03 in the paramagnetic state which pro- 
vide a basis for interpretation of 139 La NMR measure- 
ments in the manganites. The spectrum is a stan- 
dard powder pattern broadened to a significant degree 
by static lattice distortions and to a lesser degree by 
a distribution of susceptibility among lanthanum sites. 
The static lattice distortions may be due to defects or 
inherent disorder due to Jahn- Teller distortions or cal- 
cium substitution. In the temperature range studied, 
static lattice distortions around the lanthanum site do 
not change. Also, changes in the hopping rate of localized 
field inhomogcncities, such as small polarons, do not af- 
fect the spectrum. Measurements of nuclear spin-lattice 
relaxation at different points across the spectrum indi- 
cate that the relaxation mechanism is magnetic, implying 



that spin-lattice relaxation probes electron spin dynam- 
ics rather than lattice dynamics. The field dependence 
of the spin-lattice relaxation reveals several unusual fea- 
tures about the magnetic fluctuations seen at the lan- 
thanum site. First, the fluctuations are extremely slow 
at all temperatures examined here. Second, the field de- 
pendence of the fluctuations crosses over from Ti oc Hq 
near the Curie temperature to Ti oc Hq at 375K. The 
physical origin of the unusual temperature and field de- 
pendence of the spin-lattice relaxation remains unclear. 
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APPENDIX: DISTINGUISHING MAGNETIC 

AND ELECTRIC FIELD GRADIENT 
MECHANISMS OF NUCLEAR RELAXATION 

In an applied magnetic field, Hq, a spin-7/2 
nucleus with magnetic quantum number m 

7.2. —5/2, . . . , 5/2, 7/2 adopts energy levels with ener- 
gies: 

E m = mj n H (Al) 

In thermal equilibrium, the population, N m , of nuclei 
with energy E m obeys the Boltzmann distribution: 

N m oc exp (A2) 

For the measurements described here, perturbations to 
the energies, E m , due to magnetic or quadrupolar shifts 
are too small to have a significant effect on the popula- 
tions of the energy levels. The quadrupolar shifts are, 
however, large enough allow us to distinguish the po- 
sition of the central transition from the satellites in our 
spectrum measurements. Our measurements actually de- 
termine the population difference between adjacent lev- 
els. For example, the size of the signal of the central 
transition (m=+l/2 to -1/2 transtion) is proportional to 
N- 1/2 -N 1/2 . 

When measuring the spin-lattice relaxation time, r, 
one disturbs the populations of the energy levels and 
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measures the rate at which the populations return to 
equilibrium. In the so-called inversion-recovery measure- 
ments described in this paper, we use radiofrequency 
pulses to exchange the populations of two adjacent en- 
ergy levels and then measure the populations of these 
levels during their return to equilibrium. Because the 
broadening of the spectrum due to the quadrupole inter- 
action is much greater than the bandwidth of the pulses, 
the chances of affecting more than one energy level of the 
same nucleus are small, even when effects due to mag- 
netic broadening or a non zero asymmetry parameter are 
taken into account. 

In the case in which relaxation occurs due to magnetic 
fluctuations, the populations of the energy levels recover 
from a disturbance according to a set of rate equations 
1W: 



dN„ 



= -W m N m + W m+1 N m+1 + W m ^N, 



(A3) 



Where W m is the sum of matrix elements of the Hamilto- 
nian describing the coupling between the magnetic fluc- 
tuations and the nucleus. Because magnetic fluctuations 
can only induce transitions with Am = ±1, W m only in- 
cludes matrix elements describing the coupling between 
levels m and mil. Defining the population difference 
between levels m and m+1: 



L m, , m-\-l 



N, 



m+l 



N„ 



(A4) 



and solving the rate equations for the central transition 
in an inversion-recovery experiment, one finds: 



over the contributions of each satellite at any particu- 
lar frequency, giving a relaxation time,T 1 a " era9e (i/) which 
depends of the frequency position on the spectrum: 



1 — 2 exp 



= l-2j> 



7l ( l /)exp(-/ractT 1 m ) 



(A7) 



where the constants c m {v) give the intensity of the tran- 
sition between levels m and m+l at the frequency v 1 a 
quantity extracted from the fit to the spectrum. 

In the case in which relaxation occurs due to electric 
field gradient fluctuations, both Am = ±1 and Am = ±2 
transtions are allowed, giving rate equations of the form: 



dN„ 
dt 



-W m N m + W m+1 N m+1 + W m -xN m _x 



+W m+2 N m+2 + W„ 



771-2-^771-2) 



(A8) 



In a typical system, the matrix elements with Am = ±1 
and Am — ±2 are nearly the same [2l| . In this case, 
the relaxation again adopts a nearly single-exponential 
form. In contrast with the case of magnetic relaxation, 
the relaxation time of the central transition is longer than 
those of the satellites. It is possible for the opposite trend 
to occur in the situation in which electric field gradient 
fluctuations dominate, but only if the Am = ±2 matrix 
elements are over an order of magnitude larger than the 
Am = ±1 matrix elements. 



"-1/2,1/2^) OC 1 



1225 
"858 



exp 



22 



exp 



56A 75 
~~) ~ 182 



12A 1 



exp 
2t 

T 



30* 

T 



(A5) 



with different coefficients and time constants for each of 
the satellites. It turns out that the signal, n mjm +i(t), 
is well-approximated by single-exponential recovery for 
each of the transitions: 



l(t) oc 1 



t 

2 exp 

1 \ r£Ta 



(A6) 



For the inversion-recovery experiment, these effective val- 
ues relaxation time, T" 1 , adopt the following values: 

m T™ 
1/2 LOO 
3/2 1.10 
5/2 1.39 
7/2 2.58 

1/2 

where we have set T 1 / to 1.00 arbitary units. The time 
constant T™ corresponding to the central transition is 
shorter than those of the satellites, and T™ = Tf m . To 
take into account the overlap of satellites, we average 
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FIG. 1. Comparison of 575K La NMR spectra measured 
at 8.19 (filled squares) and 4.00 Tesla (open circles). Lines are 
guides to the eye. The frequency axes have been shifted for 
the sake of comparison. Inset: Comparison of narrow part of 
the spectrum. The frequency scale has been multiplied by the 
applied field to show that the width of the narrow part scales 
inversely with applied field. 

FIG. 2. Comparison of calculated powder pattern with no 
adjustable parameters and the 575K, 8.19 T spectrum. Inset 
on left is a magnified view of the narrow central feature. Inset 
on the right is a magnified view of the right side of the broad 
feature. 

FIG. 3. Same as figure ^, except the calculated powder 
pattern has been broadened by a distribution of magnetic 
shifts 

FIG. 4. Same as figure ^, except the calculated powder 
pattern has been broadened by a distribution of electric field 
gradients and magnetic shifts. 

FIG. 5. Spectra measured at several temperatures. The 
areas have been normalized to account for 1/T Curie law be- 
havior of the nuclear polarization and for changes in the Q of 
the detection circuit. 

FIG. 6. Plot of (resonant frequency) /(shift) versus temper- 
ature. The line is a fit to Curie- Weiss behavior with a Curie 
temperature of 269±19K. 

FIG. 7. Room temperature spectrum (filled squares) and 
broad part of 575K spectrum (open circles) measured at 8.19 
Tesla. Frequencies have been shifted and amplitudes scaled 
for the sake of comparison. 

FIG. 8. Temperature dependence the of inverse of the 
width of the convolution gaussian used to fit the spectra. The 
line is a fit to Curie- Weiss behavior with a Curie temperature 
of 280±30K. 



FIG. 9. Fast Fourier Transforms of spin echoes measured a 
time, T ev after selectively inverting portions of the spectrum. 
Measurements shown were performed at room temperature in 
a 4 Tesla field. 

FIG. 10. Spin-lattice relaxation times measured at 575K 
at 9.40 Tesla at different positions across the spectrum. Zero 
shift corresponds to the apex of the central transition. Times 
have been normalized to the value measured at zero shift. The 
line is a fit which assumes magnetic relaxation and overlap 
among transitions as determined from the fit in figure 0. 



FIG. 11. Spin lattice relaxation time at room temperature 
versus Hq . The fit corresponds to behavior expected for a sys- 
tem with an exponential electron spin-spin correlation func- 
tion as described in the text. 
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